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ABSTRACT 

This  paper  considers  admissibility  criteria  for  non-linear  conservation 
laws  based  on  (i)  viscosity  and  (ii)  capillarity  and  viscosity,  it  is  shown 
by  means  of  specific  examples  that  while  (ii)  yields  results  consistent  with 
experiment  for  materials  exhibiting  phase  transitions/  e.g.  a  van  der  Waals 
fluid,  (i)  does  not. 
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SIGNIFICANCE  AND  EXPLANATION 


In  the  study  of  phase  transitions  in  fluids  a  typical  model  is  provided 
by  the  equilibrium  configuration  of  a  van  der  Waals  fluid.  Such  equilibria 
show  a  fluid  may  exist  in  two  phases ,  liquid  and  vapor.  In  this  paper  we 
consider  criteria  which  should  admit  the  "physically  relevant"  liquid-vapor 
shock  wave  interfaces.  It  is  shown  that  a  criterion  based  on  both  interfacial 
capillarity  and  viscosity  yields  results  consistent  with  experiment  for 
isothermal  motions. 
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THE  VISCOSITY-CAPILLARITY  ADMISSIBILITY  CRITERION 
FOR  SHOCKS  AND  PHASE  TRANSITIONS 


*  * 

R.  Hagan  and  M.  Slemrod 

To  Jerry  Ericksen  on  the  occasion  of  his  60th  birthday. 

0 .  Introduction 

In  a  recent  paper  Serrin  [1]  reconsidered  Korteweg's  theory  of 
capillarity  [2] ,  [3]  and  applied  it  to  find  conditions  for  equilibrium  of 
liquid  and  vapor  phases  in  a  van  der  Waals  fluid.  In  subsequent  papers 
Slemrod  [4] ,  [5]  extended  Serrin* s  approach  and  proposed  Korteweg's  theory  as 
a  natural  way  of  choosing  physically  meaningful  solutions  for  dynamic  changes 
of  phase  in  a  van  der  Waals  fluid,  in  turn  James  [6]  has  shown  that  the 
Korteweg  theory  is  also  a  rational  way  of  studying  dynamic  changes  of  phase  in 
non-elasticity  as  well,  e.g.  in  materials  exhibiting  "martensitic”  or  "shape 
memory"  phase  transitions  (see  Ericksen  (7]).  In  this  paper  we  carry  this 
program  one  step  further:  We  propose  the  Korteweg  theory  to  provide  the 
universal  admissibility  criterion  for  isothermal  motions  of  compressible 
thermo-elastic  fluids.  Indeed  it  is  our  belief  that  the  Korteweg  theory 
combined  with  the  introduction  of  the  natural  thermal  dissipation  of  heat 
conduction  should  provide  the  universal  admissibility  criterion  for 
compressible  thermo-elastic  fluids  in  general.  For  the  purposes  of  this 
paper,  however,  we  restrict  ourselves  to  the  isothermal  case.  We  shall  show 
by  specific  examples  when  the  theory  agrees  and  disagrees  with  predictions 
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based  on  viscosity  alone.  Furthermore  we  will  point  out  the  power  of  the 
Korteweg  theory  in  choosing  admissible  singular  surfaces  for  conservation  laws 
which  are  neither  genuinely  non-linear  nor  even  hyperbolic,  e.  g.  motion  of  a 
van  der  waals  fluid. 

The  paper  is  divided  into  four  sections.  The  first  section  derives  the 
balance  laws  of  compressible  fluid  mechanics  in  Lagrangian  form.  The  second 
section  discusses  admissibility  with  respect  to  viscosity  alone.  The  third 
section  considers  admissibility  with  respect  to  the  viscosity-capillarity 
condition  as  derived  from  Korteweg* s  theory.  In  particular  we  show  that  the 
usually  accepted  ideas  on  admissibility  of  shocks  in  ideal  gases  are 
consistent  with  the  viscosity-capillarity  condition.  Finally  the  fourth 
section  applies  the  viscosity-capillarity  condition  the  the  issue  of  phase 
transitions  in  a  van  der  Waals  fluid. 

We  note  that  recently  Shearer  [19]  has  also  made  a  study  of  the  viscosity 
and  viscosity-capillarity  criteria.  Basically  his  work  centers  on  what  is 
here  given  in  Examples  5  and  6  of  Section  3.  Aside  from  this  overlap  our 
paper  is  devoted  to  issues  not  considered  in  [19]. 


1.  Preliminary  discussion  of  the  equations  of  compressible  fluid  flow. 

We  follow  the  presentation  of  Courant  and  Friedrichs  [8]  of  a  Lagrangian 
description  based  on  the  law  of  conservation  of  mass.  The  fluid  flow  is 
thought  of  as  taking  place  in  a  tube  of  unit  cross  section  along  the  x-axis. 

We  attach  the  value  X  ■  0  to  any  definite  "zero"  section  moving  with  the 
fluid.  For  any  other  section  we  let  X  be  equal  in  magnitude  to  the  mass  of 
the  fluid  in  the  tube  of  unit  cross  sectional  area  between  that  section  and 
the  zero  section.  Analytically  the  quantity  X  satisfies  the  relation 

x "  ^((o!t)  p(x't)dx  •  <*-i> 

Here  p(x,t)  denotes  the  density  of  position  x  and  time  t  and 
x  »  x(X,t)  denotes  the  position  of  the  particle  so  that  a  mass  X  of  fluid 
is  enclosed  in  the  tube  bounded  by  x(X«t)  and  x(0,t).  Differentiation  of 
(1.1)  implies  1  -  Xx<X,t)p (x <X,t) ,t) .  Let  p(x(X,t),t)  -P(X,t),  w(X,t)  « 

*v  .1  • 

p(X,t)  (the  specific  volume),  x(x,t)  *  u(X,t)  *  Xt(X,t)  (the  velocity). 

We  denote  the  stress  by  t,  the  specific  internal  energy  by  e,  the 
specific  heat  absorption  by  q,  the  heat  flux  by  h,  the  specific  body  force 
by  b.  Then  the  equations  of  balance  of  linear  momentum,  energy,  and  mass 
become 


PX  »  +  Pb  , 

pc  =  Pq  +  TXx  +  hx  , 


p  +  pxx  ■  0  , 


d 


where  •  «  .  We  now  apply  the  chain  rule  and  rewrite  this  system  in  the 

Gt 


terms  of  the  independent  variables  X,  t  to  obtain 


*tt  -  Tx  *  b  ' 


e  *  q  +  ty  +  h 

£  'l  A  4-V  X 


ltx  T  nx  ' 


(1.2) 


(PXX>t  = 


0 


where  we  have  used  the  fact  that  P(x,t)  ■  Pt(X#t)«  The  third  of  these 
equations  is  automatically  satisfied  since  pxx  *  1. 

The  above  set  of  balance  laws  must  be  supplemented  by  constitutive 
relations.  We  assume  the  fluid  is  compressible,  isotropic,  and  thermo-elastic 
where  the  stress,  internal  energy,  and  heat  flux  satisfy 

T  =  -W (w,T)  , 


e  -  e(w,T)  , 
h  -  h(w,T)  , 
S  *  s(w,T)  , 


(1.3) 


*  -  i|i(w,T)  , 

where  *  (the  pressure)  and  s  (the  specific  entropy)  are  deriveable  from 


the  Helmholtz  free  energy  potential 


a*  a* 

"  “  ai  '  8  "  "fr  * 


(1.4) 


We  now  make  the  simplifying  hypothesis  that  the  fluid  is  imbedded  in  a 
"heat  bath"  so  that  the  motion  is  Isothermal  (T  -  positive  constant)  and  that 
there  are  no  body  forces.  Mathematically  this  means  that  q  is  assumed  to  be 
adjusted  so  that  (1.2a)  is  always  satisfied  identically  with  T  -  constant  and 
b  5  0. 

If  the  motion  is  isothermal  * (w,T)  is  a  function  of  w  alone  and  we 
set  *(w,T)  »  p(w).  Hence  (1.2)  is  equivalent  to  the  first  order  system 


ufc  "  -p(w)x  * 


w  *  u  « 
t  x 


(1.5) 


where  we  assume  p  €  C  (0,°*).  As  is  easily  seen  (1.5)  is  either  hyperbolic  or 
elliptic  depending  on  the  sign  of  p’(w). 

A  curve  F  :  X  =  o(t)  across  which  u,  w  experience  jumps  is 

called  a  singular  surface  or  shock  wave.  If  T  is  a  singular  surface  let 


(0<t),  t)  be  a  fixed  point  on  the  graph  of  T  and  let  U  =  o( t).  Denote 
by  u+,  w+,  u_,  w_  the  respective  limits  from  the  right  and  the  left  as 
(X,t)  ♦  (0(1),  t)  for  (u,w).  If  we  put  [uj  *  u+  -  u_,  [w]  =  w+  -  w_#  < 
then  we  known  the  Rankine-Huginot  jump  conditions  must  be  satisfied  across 
r,  i.e. 


U[u]  -  [p]  *  0 
U[w]  +  [u]  ■  0 


( 


o 

Of  course  (1.6)  implies  U  satisfies  the  equation 


2.  Admissibility  with  respect  to  the  viscosity  criterion 


Let  us  imbed  the  inviscid  equations  (1.5)  in  a  viscous  formulation,  i.e. 

we  take 


T  -  -p(w)  +  y*^  (2.1) 

where  y  >  0  is  the  (assumed  constant)  viscosity. 

Now  let  T  :  X  »  o( t)  be  a  singular  surface  for  (1.6).  We  ask  the 
question:  Are  solutions  of  (1.5)  in  the  neighborhood  of  the  singular  surface 
T  limits  of  solutions  ( 1 . 2a ) ,  ( 2. 1 )  as  y  +  0+?  While  this  problem  has  a 
long  history  dating  back  to  Rayleigh  [9]  it  is  the  more  recent  discussions  of 
Wendroff  [10]  and  Dafermos  [11]  we  shall  follow. 

Let  (0(t),  t)  be  a  fixed  point  on  the  graph  of  T  and  let  u+,  w+,  u_, 
u  be  as  in  Section  1.  We  look  for  a  traveling  wave  solution  of  (1.2), 
(2.1)  given  by 


u(X,t)  *  u(5),  w(X,t)  -  w(€),  5 
It  follows  that  u,  w  must  satisfy 

-Uu'  =  (-p  +  u' ) '  , 
-U$'  *  u'  , 


X-Ut 

P 


(2.2) 


where  '  *  ^  . 

In  order  for  u,  w  to  approximate  the  discontinuous  profile  of  the 
solution  to  (1.5)  we  require 

(u(-"),  w(-*»),  u(-H»),  w(+*°))  *  (u_,  w_,  u+,  w+)  .  (2.3) 

Definition  2. 1 .  If  there  exists  u,  w,  C1  functions  so  that  (2.2),  (2.3)  are 
satisfied  for  all  points  (0 ( t),  t)  we  say  the  singular  surface  satisfies  the 
viscosity  admissibility  criterion. 

Theorem  2. 2.  The  viscosity  admissibility  criterion  is  satisfied  if  and  only 
if 


,  -p(w)  +  p(w  ) 

-0  +  ( - - - — )  >  0  (<  0) 


if  U  >  0(<  0) 


(2.4) 


for  every  value  w  between  w_  and  w+. 

In  other  words  for  (w+  -  w_)U  >  0  ( (w+  -  w_)U  <  0)  the  chord  which 
joins  (w_,  p(w_)  to  (w+,  p(w+))  lies  above  (below)  the  graph  of  the 
function  p(w)  for  w  between  w_  and  w+. 

Proof.  Integrate  (2.2)  from  -*  to  £ .  It  then  easily  follows  that  w(£) 
satisfies  the  first  order  equation 

U2(w(£ )  -  w_)  +  p(w(£))  -  p(w_)  =  -Uw’(£)  •  (2.5 

If  U  /  0  and  w_  <  w+(w_  >  w+)  then  we  must  have  w'(£)  >  0  (w'(£)  <  0) 
and  hence  (2.4)  follows. 

C 

In  the  exceptional  case  of  a  static  singular  surface  where  0-0  we 
have  £  =*  X/u  and  (2.2)  becomes  p(w(£))'  "  0.  So  if  w  is  a  C1  function 
satisfying  w(-*)  =■  w,  w(-*“)  -  w+  we  see  p(w(£))  =  p(w_)  ”  p(w+)  for  all 
-*»<£<«  and  hence  p(w)  must  equal  p(w_)  =  p(w+)  for  all  w  between 
w_  and  w+.  So  unless  p  is  identicaly  constant  on  [w_,  w+]  a  static 
singular  surface  is  never  admissible  according  to  the  viscosity  criterion. 
E*»mple  1.  Let  p(w)  have  the  graph  shown  in  Figure  1  where  p'  <  0,  p"  >  0 
for  w  >  w+.  In  this  case  the  chord  connecting  (w„,  p(w_))  to  (w+/  p(w+)) 
lies  above  the  graph  of  p(w)  and  the  singular  surface  T  :  X  =  Ut  with 

/-p(w  )  +  p(w_) 

- —  ( 2.6 

w+-w 


is  admissible  according  to  the  viscosity  criterion. 


Figure  1 


Example  2.  Let  p(w)  have  the  graph  shown  in  Figure  2.  In  Figure  2  the 
chord  joining  (w_,  p(w_ ) )  to  (w+/  p(w+))  again  lies  above  the  graph  of 
p(w)  for  w_  <  w  <  w+  even  though  the  chord  is  tangent  to  the  graph  at 
(w+,  p(w+)).  Thus  the  singular  surface  T  :  X  *  Ut  with  u  given  by  (2.6) 
is  again  admissible  according  to  the  viscosity  criterion.  The  purpose  of  this 
example  is  to  point  out  that  the  viscosity  criterion  allows  for  tangency 


constructions . 


3.  The  viscosity-capillarity  criterion 


P 


i 

i 


'i 


4 


4 


The  problem  with  the  viscosity  criterion  is  that  it  allows  only  the 
viscosity  as  a  higher  order  singular  correction  in  the  thermo-elastic  fluid's 
constitutive  relation.  In  line  with  the  generally  held  belief  that 
consideration  of  other  physically  meaningful  effects  should  help  pick  out  the 
physically  relevant  solutions.  Slemrod  [4]  suggested  a  new  condition  termed 
the  viscosity-capillarity  criterion  based  on  Korteweg's  theory  of  capillarity 
([2],  [3]).  In  this  section  we  will  review  the  viscosity-capillarity 
criterion  and  apply  it  to  some  examples. 

According  to  Korteweg's  theory  the  stress  T  in  our  one-dimensional, 
isotropic,  isothermal  formulation  will  be  given  by 

x  -  p(w)  +  B(w)w2  -  Clwjw^  +  W(w)ux  .  (3.1) 

For  simplicity  we  shall  consider  only  the  constant  coefficient  case 

U(w)  =  y  ,  a  positive  constant, 

2 

C(w)  »  y  A  ,  A  a  constant, 

B(w)  =  -y  ,  D  a  constant. 

In  this  special  case  (1.2a),  (3.1)  become 

2  2  2 

ut  =  (-p(w)  -  y  Aw^  +  y^  -  y  ^V'x  ' 

(3.2) 

wt  =  "X  ‘ 

Now  again  let  T  be  a  singular  surface  with  o(t),  U,  w_,  w+,  u_,  u+  as 
given  in  Section  1.  Again  we  look  for  traveling  wave  solutions 
u(X,t)  *  u(5),  w(X,t)  =  w(5),  £  =  , 

this  time  for  equations  (3.2).  It  then  follows  that  u,  w  must  satisfy 

-Uu*  *  (-p(w)  -  Aw"  -  Dw'2  +  u')'  , 

A  A  (3*3) 

-Uw'  =  u'  . 


I 
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In  order  that  u,  0  will  approximate  the  discontinuous  profile  we 
require  again  that 

(u(-"),  w(-*>),  u(4-»),  w(+“))  =■  (u_,  w_,  u+,  w+)  .  (3.4) 

These  considerations  motivate  the  following  definition,  where  from  now  on 
we  drop  the  A  superscript. 

Definition  3.1.  If  there  exist  u  e  C2,  w  €  C*  satisfying  (3.3),  (3.4)  for 
all  points  «J(t),  t)  we  will  say  the  singular  surface  T  satisfies  the 
viscosity-capillarity  admissibility  criterion  (A,D). 

Since  the  criterion  may  possibly  be  satisfied  for  one  (A,D)  pair  and  not 
another  the  dependence  of  the  criterion  on  A  and  D  has  been  explicitly 
noted  in  the  above  definition. 

Lemma  3. 2.  The  viscosity-capillarity  criterion  (A,D)  is  satisfied  if  and  only 

if  there  exists  w (5 )  €  a  solution  of  the  second  order  equation 

Aw"  +  Uw'  +  Dw'2  +  f (w,U)  -  0  ,  (3.5) 

w(-"»)  *  w  ,  w(+»)  =  w  ,  (3.6) 

-  + 

where 

f(C,U)  -  U2(5-w_)  +  p(C )  -  p(w_)  .  (3.7) 

Proof.  Integrate  (3.3)  from  -*®  to  5. 

□ 


Theorem  3.3.  Assume  (3.5)  possesses  precisely  two  equilibrium  point  w  =  w_, 
w'  ■  0  and  w  *  w+,  w'  =  0  with  w_  <  w+,  p'(w_)  <  0,  p' (w+)  <  0.  Then 
there  is  a  solution  w(€)  of  (3.5)  connecting  these  equilibrium  points  if  the 
following  holds: 

(i)  A  >  0,  U  >  0  and 

(a)  (-p' (w+)  )^2  <  u  <  (-p‘ (w_) )^2  , 

(b)  the  chord  connecting  (w_,  p(w_))  and  (w+,  p(w+))  lies  above 
the  graph  of  p(w)  for  w_  <  w  <  w+  , 


I 


(c)  for  some  4  >  w+,  exp(—  C)fU,U)dC  ■  0  ; 

(ii)  A  >  0,  U  <  0  and 

(a)  -(-p'(w+))1/2  <  U  <  -(-p'fwj)1^  , 

(b)  the  chord  connecting  (w_#  p(w_))  to  (w+,  p(w+))  lies  below  the 
graph  of  p(w)  for  w_  <  w  <  w+  , 

w 

(c)  for  some  4  <  w_,  /&+  exp(jp  C)f(C,U)dC  *  0  t 
(Hi)  A  <  0,  V  <  0  and 

(a)  <-p'(w+))1/2  <  -u  <  (-p'tw.))1^  f 

and  i-b,  ii-c  hold; 

(iv)  A  <  0/  U  >  0  and 

(a)  -(-p* (w+) fc  <  -u  <  -(-p* (w_) )^2  , 

and  ii-b,  i-c  hold. 

Proof.  Rewrite  (3.5)  in  first  order  form 


and  set 


w*  »  v  , 

Av*  =  -f (w,U)  -  Uv  -  Dv2  , 

2D  2 

H(WfV)  =  exp(~  w)v  +  G(w,U),  where 


(3.8) 


G(w,U)  *  |  JJJ  exp(jp  C)f(C,U)dC  . 
Let  a  *  { (w, v) ;  w_  <  w  <  4,  H(w,v)  <  H(w_,0)> 


*  {  (w *  v)  j  w_  <  w  <  4,  exp(^  w)v2  <  ~  /  exp(-^  ?)f(?,U)dC 

A  AW  A 

where  4  is  now  assumed  to  be  the  smallest  such  4  so  that  i-c  holds.  From 
i-b  and  i-c  we  see  a  is  a  simply  connected  bounded  open  set  in  R2 

containing  (w+,0).  Furthermore 

dH  2  2D  .  2  .  2 

~  =  -  —  U  exp(j-  w)v  *  -const,  v  for  const.  >  0 

where  (w,v)  €  a.  Thus  any  orbit  which  enters  Q  remains  in  Q.  Also  i-a 


-12- 


shows  (w_,0)  is  a  saddle  and  an  elementary  analysis  of  this  saddle  shows 

that  it  possesses  an  unstable  manifold  which  enters  ft.  By  LaSalle's 

invariance  principle  ([12]#  p.  316)  an  orbit  entering  ft  on  this  manifold 

must  approach  an  equilibrium  point  since  H  is  a  Liapunov  function  on  ft. 

dH 

Such  an  orbit  cannot  approach  (w_,0 )  since  <  0  in  ft  and 
H(w#v)  <  H(w_,0)  for  (w,v)  in  ft#  hence  it  approaches  (w+,0). 

(ii)  The  proof  here  is  analagous  to  (i)  with 

G(w,U )  »  7 /"  exp(“  C)f(C,U)dC  . 

The  saddle  now  is  at  (w_#0)  and  ft  “  {(w,v)»  t  <  w  <  w+,  H(w#v)  <  H(w+,0)} 
where  i  is  assumed  to  be  the  largest  such  t  so  that  ii-c  holds. 

(iii)  The  proof  is  analagous  to  (i)  with  G  and  ft  as  given  in  (ii). 

(iv)  Again  the  proof  is  analagous  to  (i)  with  G  and  ft  as  given  in  (i). 


Corollary  3.4.  (i)  The  conclusion  of  Theorem  3.3(i)  holds  if  i-c  is  replaced 

by  (i-c’)s  p"  >  0  for  w  >  w+. 

(ii)  The  conclusion  of  Theorem  3.3(ii)  holds  if  ii-c  is 
replaced  by  ( ii-c ' ) :  p"  <  0  for  w  <  w_ . 

Proof,  (i)  We  note  (i-c')  implies  p'(w)  -  p' (w+)  >  0  for  w  >  w+  and  hence 


by  i-a  p' (w)  +  U2  >  e  for  some  c  >  0.  Thus  p(w)  -  p(w+)  + 

w_ 

U2(w-w.)  >  e(w-w.)  for  w  >  w .  and  /  exp(—  C)f(C#U)d£  ♦  -*  as 

T  +  T  W  A 

*•*■■*».  Hence  there  is  an  l  so  that 

w 

Jt”  exp(^  ?  )f(?,U)dC  -  0,  l  >  w+  . 

(ii)  The  proof  is  similar  to  (i). 


□ 


Corollary  3.5.  If  p'  <  0,  A  >  0  and  i-a  holds  then  the  conclusion  of 


Theorem  3.3(1)  holds 


Proof.  Apply  Corollary  3.4 


□ 

w_ 

Corollary  3.6.  If  p'  <  0,  p"  >  0,  A  <  0  and  lim  /  p(C)dC  ■  +•  then  the 

w+0+  w 

conclusion  of  Theorem  3.3(lil)  holds. 

Proof.  Here  i-b,  i-d  are  automatically  satisfied  and  the  rest  is  an  obvious 
application  of  Theorem  3.3(iii). 

□ 

Example  3.  Consider  an  ideal  gas  with  constitutive  relation  *(w,T)  *  RTw  1 
where  R  is  a  positive  constant.  As  the  chord  joining  (w_,p(w_) )  to 
(w+/p(w+) )  lies  above  the  graph  of  p(w)  between  w_  and  w+,  the  singular 
surface  T  ;  X  ■  Ut  with  U  given  by  (2.6)  is  admissible  according  to  the 
viscosity  condition.  Corollary  3.5  shows  T  is  also  admissible  according  to 
the  viscosity-capillarity  condition  (A,D)  for  A  >  0.  On  the  other  hand 


Oorollary  3.6  shows  that  T 


is  admissible  according  to  the  viscosity-capillarity  condition  (A,D)  if 
A  <  0.  Thus  A  <  0  admits  rarefaction  shocks  in  an  ideal  gas.  For  this 
reason  we  view  A  >  0  as  the  only  physically  reasonable  choice. 

Example  4.  Let  p(w)  be  as  in  Figure  1.  In  Example  1  we  saw  T  :  X  *  Ut 
was  admissible  according  to  the  viscosity  criterion  where  U  is  given  by 
(2.7).  Now  we  see  by  Corollary  3.4 (i)  that  this  same  T  is  also  admissible 
according  to  the  viscosity-capillarity  criterion  (A,D),  A  >  0. 

Example  5.  Assume  p(w)  has  the  graph  as  shown  in  Figure  2.  In  this  case 

/w 

f(C,U)d?  has  the  graph  shown  in  Figure 

w+ 

3. 
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Figure  3 


We  note  the  point  (w+,  G(w+))  is  a  point  of  inflection.  Typical  phase 
flows  exiting  from  w  -  w_,  v  *  0  will  generally  be  unbounded  and  only  in 
exceptional  circumstances  will  we  have  w(+“»)  -  w+.  of  course  (3.5),  (3.6) 
has  the  rather  simple  mechanical  analogue.  Namely  one  visualizes  a  bead 
rolling  down  the  hill  given  by  the  graph  of  G(w).  Admissibility  according  to 
the  viscosity-capillarity  criterion  (A  >  0,  D  *  0)  would  require  the  damping 
force  of  rolling  motion  given  by  Uw'  sufficient  to  bring  the  bead  to  rest  at 
the  inflection  point  (w+,  G(w+))  -  a  highly  unlikely  circumstance.  Hence 
except  for  this  exceptional  case  the  singular  surface  T  :  X  *  Ut  where  U 
is  given  by  (2.6)  will  be  admissible  according  to  the  viscosity  criterion  and 
inadmissible  according  to  the  viscosity-capillarity  criterion  (A  >  0,  D  *  0). 
Example  6.  Let  p(w)  be  as  shown  in  Figure  5. 


Figure  4 


In  Figure  4  the  chord  joining  (w_,  p(w_) )  to  (w+,  p(w+))  again  lies 
above  the  graph  of  p(w)  for  w_  <  w  <  w+  even  though  the  chord  is  tangent 
to  the  graph  at  (w+l  p(w+) ) .  So  the  singular  surface  T  s  X  ■  Ot  with  U 
given  by  ( 2.7)  is  admissible  according  to  the  viscosity  criterion*  However  if 
we  plot  the  potential  G(w)  *  /*  f(C,U)dC  we  find  G(w)  has  the  graph  shown 


in  Figure  5 


Here  again  as  in  Exanqple  4  (w+,  G(w+))  is  an  inflection  point.  Now 
however  instead  of  the  typical  flows  of  (3.5)  exiting  w  =  w_,  w'  -  0  being 
unbounded  these  flows  will  either  be  unbounded  or  approach  the  equilibrium 
point  w  ■  v*,  w'  *  0,  modulo  the  exceptional  case  that  the  f(w) 
approaches  w  *  w*,  w*  -  0.  So  again  the  singular  surface  r  ;  x  «  Ut  with 
U  given  by  (2.6)  is  admissible  according  to  the  viscosity  criterion  and 
inadmissible  (modulo  an  exceptional  case)  according  to  the  viscosity 
capillarity  condition  (A  >  0,  D  *  0). 


4*  Dynamics  of  phase  transitions  in  a  van  der  Waals  fluid 


In  this  section  we  examine  the  admissibility  of  singular  surfaces  for 
isothermal  motions  of  a  van  der  Waals  fluid.  In  a  van  der  Waals  fluid  the 
constitutive  relation  for  the  pressure  is  given  by 

OT  m 

n  (w,T)  »  ~—r - r  ,  0  <  b  <  w  <  • 

W-D  2 
W 

where  a,  b,  R  are  positive  constants  (see  for  example  [13] f  [14],  [15]). 
Again  we  set  p(w)  -  H(w,t)  for  isothermal  motions,  T  »  positive  constant. 
For  T  sufficiently  small  p(w)  has  a  graph  as  shown  in  Figure  6. 


Figure  6 

Actually  in  what  follows  we  shall  not  require  anything  so  specific  as  the 
van  der  Waals  equation  of  state.  Instead  we  assume  p(w)  satisfies  the 
following  hypotheses: 

(Hi)  p' (w)  <  0  for  w  €  (b,a)  (@,«),  b  >  0  , 

(H2)  p'  (a )  -  p'tS)  -  0  , 


(H3)  p'  (w)  >  0  for  w  €  (0,8)  where  b<a<8<co. 


(H4 )  p(w)  >  0  for  w  €  (b,*0)  , 
.2, 


(H5)  p  €  C  (b,°» )  , 

(H6)  p(w)  <  p(o)  as  w  ♦  •  , 

(H7)  p(w)  has  one  inflection  point  in  (8,*)  , 

(H8)  p(Y )  -  p(B),  p(5)  -  p<a),  b<Y<a,  8<«<-. 

The  domains  (b,a)  and  (8#“)  will  be  called  the  a  (liquid)  phase  and 
8  (vapor)  phase  respectively. 

We  shall  now  consider  the  admissibility  of  several  singular  surfaces 
separating  liquid  and  vapor  phases  of  a  van  der  Waals  fluid. 

Example  7.  When  w^  and  w^  are  as  shown  in  Figure  7  the  chord 
connecting  (Wj  ,p(w1  ))  and  (Wj  ,p(w^  ))  lies  above  the  graph  of  p(w) 


between  w^  and  w^  as  so  r  :  X  »  O^t  with 


/-p(w1  )+p(w1  ) 
- 1 - — 


w  -w 
1  1 
+ 


is  admissible  according  to  the  viscosity  criterion.  Theorem  3.3(i)  shows  T 
is  also  admissible  according  to  the  viscosity-capillarity  condition  (A,D), 


A  >  0. 
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Figure 


Example  8.  When  w2  and  w2  are  shown  in  Figure  7  the  same  chord-graph 
argument  as  in  Example  7  shows  F  :  X  *  U^t  with 


/^r 

u  -  +  / - 


)+p(»2  ) 


- 

2+  2- 

is  admissible  according  to  the  viscosity  criterion.  However  an  argument 
analagous  the  one  of  Example  6  shows  that  (modulo  an  exceptional  circumstance) 
T  will  not  be  admissible  according  to  the  viscosity-capillarity  condition 
(A,D),  A  >  0,  D  *  0.  We  expect  that  admissibility  according  to  the  viscosity- 
capillarity  criterion  will  not  be  aided  if  0^0. 

Example  9.  For  and  w3  shown  in  Figure  7  the  remarks  made  in  Example 

8  again  apply.  So  T  :  X  “  Ujt  with 


/-P<V>*P(»3 
o -  -  / - - - : 


3+  ’3. 


is  admissible  according  to  the  viscosity  criterion  but  not  generally 
admissible  according  to  the  viscosity-capillarity  condition  (A,D),A  >  0. 
Example  10.  For  w4  and  w4+  shown  in  Figure  7  the  singular  surface  is 
F  :  X  ■  0  with  U4  «  0  since  p(w4  )  “  p(w4  ).  By  the  remark  in  Section  2 
T  cannot  be  admissible  according  to  the  viscosity  criterion.  On  the  other 
hand  an  elementary  quadrature  solution  of  (3.5),  (3.6)  shows  that  if  w4 
and  w4  satisfy 


Jw  +  exp(+  C)(p(C)  -  p(v4  ))<K  -  0 


-2 


then  T  is  admissible  according  to  the  viscosity-capillarity  condition 
(A,D) •  Of  course  if  D  -  0  the  above  relation  reduces  to  the  well  known 
Maxwell  equal  area  rule  [1]/  [13],  [14].  If  D  0  the  relation  falls  within 
the  new  class  of  "rules"  proposed  by  Serrin  [ 1 ] . 

We  note  that  for  non-zero  U  when  the  chord  connecting  (w_,p(w„)) 
and  (w,,p(w,))  cuts  the  graph  of  p(w)  at  intermediate  value  of  w, 

I 

Theorem  2.2  precludes  any  such  construction  being  admissible  according  to  the 
viscosity  criterion.  Such  a  construction  may  however  be  admissible  according 
to  the  viscosity-capillarity  condition  (A,D) .  We  now  examine  this  situation 
in  detail . 

example  11.  Let  w_  be  given  (b,a)  and  w+(U)  be  a  solution  in  (0,*°)  of 
^  the  Rankine-Huginot  relation  with  U  >  0,  i.e.  w+(U)  satisfies 

/-p(w+(U) )+p(w  ) 

°-+/ - ^(b)-w_  '  * 

I  The  line  T  :  X  =  Ut  defines  a  singular  surface.  We  consider  the 

admissibility  of  T  according  to  the  viscosity-capillarity  condition  (A,D) 

2 

when  f(?,U)  =  p(C)  -  p(w_)  +  U  (£-w_)  possesses  four  zeroes,  namely  w_, 

O  ♦ 

w+(U)  of  course,  plus  two  additional  zeroes  w(U),  w  (U)  where 

O  *  h 

o  <  w(U)  <  0  <  w  (U)  <  w+(U).  Thus  w  (U)  is  the  smallest  solution  of 

f(C»U)  *  0  in  (0,00)  and  w+(U)  is  the  solution  of  f(C,U)  =0  in  (0,®) 

2  2 

which  satisfies  p'(C)  +  U  >0.  We  also  assume  U  <  -p'(w_)  (see  Figure 

8). 
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b 


a 


w(U)0  w*(U)  6 


w+(U)\ 


Figure  8 


Again  the  admissibility  of  T  according  to  the  viscosity-capillarity 
criterion  (A,D)  is  equivalent  to  the  existence  of  a  solution  of  (3.5),  (3.6) 
or  equivalently  (3.8)  with  boundary  conditions  (w(+°“),v(-m>»)  )  *  (w+(U),0), 
(w(-“),v(-«)  =  (w_,0).  We  now  state  our  main  theorem. 

Theorem  4.1.  (i)  Let  A  >  Then  there  exists  U  >  0  so  that  the  singular 

•  11  O 

surface  T  j  X  =  Ut  for  the  solution 


w  =  w  (U) 

X  <  Ut  ,  X  >  Ut 

u  =  u  -  U(w  (U)  -  w  ) 


u 


(4.1 


of  (1.5)  is  admissible  according  to  the  viscosity-capillarity  condition  (A,D) 


if  either 

* 

(I)  a  >  w  >  Y  and  /W  (0)  exp(-^  ?  )f  (?  ,0 )dC  <  0 
or 

(I*)  b  <  w_  <  Y  , 

and 

(II)  there  exists  UM  >  0  so  that 

* 

w  (U  ) 

Jw  M  exp(|p  C)f(C,UM)dC  «  0  . 

(ii)  Let  A  >  ~.  Assume  (I)  or  (I*)  and  (II)  of  (i)  above  hold.  Then  for 
8 

all  U,  0  <  U  <  U,  the  singular  surface  T  s  X  ■  Ut  for  the  solution 


(w  «  w  w  «  w  (U) 

X  <  Ot  ,  X  >  UT  (4.2) 

u  *  u_  u  **  u_  -U(w  (U)  -  w_ ) 


of  (1.5)  is  admissible  according  to  viscosity-capillarity  criterion  C£.,D),  U 
is  as  given  in  (i). 

The  theorem  will  be  proven  via  a  sequence  of  lemmas. 

Lemma  4. 2.  Let  U  >  0,  A  >  0.  Then  the  equilibrium  point  (w_,0)  of  (3.8) 
is  a  saddle.  Futhermore  Yy(w_,0),  the  orbit  exiting  from  the  saddle  in 
the  v  >  0  half  plane,  either  approaches  an  equilibrium  point  of  crosses 
the  w_  axis  for  increasing  5 • 

Proof.  Set  n (w,U)  =  w' (5).  From  (3.5)  we  find  n  satisfies 

A  n  £  +  un  +  Dn2  +  p(w)  -  p(w_)  +  u2(w-w_)  =  0  ,  (4.3) 

n<w_,0)  -  0  ,  (4.4) 

and  hence 


-24- 


(4.5) 


n (w,(J)  -  -  I  u  Jw  exp(-^  (?-w))n(?,o)d? 

AW  A 

-  I /w  exp(^  (C-w) )f (S ,U)d£  . 

AW  A 


From  (Hi)  -  (H6)  we  know 

f(C,U)  >  0  on  (b,w_)  u  (w(U),w  (0))  u  (w+(U),-)  , 

<  0  on  (w_,w(U))  u  (w  (U),w+(U)) 


If  the  solution  exiting  from  the  saddle  doesn't  cross  the  w  axis  or  approach 
an  equilibrium  point  it  must  tend  to  +•  as  £♦+•.  But  for  such  a 
solution  we  know  from  (4.5)  that 

n2(w,U)  <  -  e*p(jp  (C-w))f(C,U)dC  . 

So  by  L'Hopital's  rule  lim  h2(w,U)  <  -  ~  lim  f(w,U)  **  -  «*.  Hence  w 

*r*+  <•  W+-+  - 

cannot  approach  +  “  and  the  lemma  is  proven. 


□ 


+  ~ 

Definition  4.3.  If  Yu(w_,0)  intersects  the  w-axis  define  w(0)  to  be  that 

value  of  w  so  that  (w(U),0)  e  ljlw_,0)  and  w(U)  is  the  first  point  of 

crossing  the  w-axis  of  Y*(w_,0).  If  1y(w_,0)  doesn't  intersect  the  w  axis 

+ 

define  (w(U),0)  to  be  the  equilibrium  point  approached  by  Yyfw^O)  as 
£  ♦  +  *. 


Lemma  4.4.  Let  n(w,U)  be  the  solution  of  (4.3),  (4.4).  Then  for  A  >  — , 

0  <  U|  <  U2,  we  have  h(w,U.j)  >  h((w,U2)  on  (w_,w(U2)).  Furthermore  we  have 
w(U2)  <  wU^). 

Proof.  A  direct  calculation  from  (4.3),  (4.4)  shows 


3  a 


i 


Sw  ntw-,D)  --S* 


( 1-4A )U 


A ( U ^-4A ( p ’  (w_)+U2))1/2 


2 

Since  p'(w_)  +  U  <  0  we  see 
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|5|:n(«.,0)  <  4  if  1-4A  >  0 


<  -  —  if  1-4A  <  0  . 

2A 

Hence  |^j  n(w_,U)  <0  if  A  >  “> 

«v  <v 

Now  suppose  iKw/U^)  -  h(w,u2)  has  a  zero  on  (w_,w(U^))  n  (w_,w(u2>). 
Let  v  be  the  smallest  such  zero.  Then  from  (4.3)  we  have 

A  (v,^)  -  (v,u2»  -  n(v,u1)(o2-o1)  +  (v-wj (u2-u2)  . 

Since  U2  >  U1  >  0,  v  >  w_,  and  n(v,U1)  >  0,  we  have 

A  (v.O,)  -  (vfU1)J  >  0.  Hence  ~  OM^)  >  ^  (v,U2)  which 

contradicts  the  assumption  that  v  is  the  smallest  zero  of 
T)(w(U1)  -  n(w,U2).  Thus  ntWfO^  >  n<w,u2)  on  (w_,w(U1>)  n  (w_,w(U2>)  and 
hence  w(U2>  <  w(U^)« 

Definition  4.5.  For  the  first  order  system  (3.8)  set 
G(w,U)  (u)  exp(^  C)f(C,0)dC  , 

H(w,v)  -  exp(~  w)v2  +  G(w,U),  and 

fly  -  { (w,v) ;  w*(U)  <  w,  H(w,v)  <  H(w*(U),0)}  . 

Lemma  4.6.  fl  is  a  simply  connected  bounded  open  set.  Furthermore  if 

\ 

A  >  0,  0  >  0,  w* (U)  <  w+(U)  any  orbit  of  (3.8)  entering  fly  has  (w+(U),0) 
for  an  w-limit  set. 

Proof.  From  (H3),  (H4),  and  (H7)  we  see  there  is  a  unique  i  >  w*(U)  so  that 

/  ,  exp{~  C )f (C ,U)cK  ■  0.  Hence  A„  is  simply  connected  bounded  and 

w*(U)  A  _  u 

,  4  <JH  2v  „  .  2D  ,  dH  ,  2  _ 

open.  Also  since  -gr  * - —  U  exp(—  w)  we  see  ^  -  const,  v  for 

const.  >  0  inside  fly. 

Because  any  orbit  entering  fly  remains  in  fly  LaSalle's  invariance 
principle  [12)  asserts  that  orbit  approaches  the  only  equilibrium  point  in 


fly,  namely  (w+(U),0), 


1 
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Lemma  4.7.  Let  A  >  0/  U  >  0,  w*(U)  <  w+(U).  Then  the  solution  of  (3.8) 
exiting  from  the  saddle  (w*(U),0)  in  the  v  >  0  half  plane  enters  8^  and 
has  (w+(U),0)  for  an  w-limit  set. 

Proof .  Linearize  about  (w*(u),0)  and  use  Lemma  4.6. 

□ 


Lemma  4.8.  For  0  <  U1  <  U2  we  have  8  c  8  . 

U2  ui 

Proof .  Note  8y  =  {(w,v);  wMUj)  <  w,  v2  <  -  j  f"*(U  )exP^  OfUrU^dC)  , 
Q0  -  t(w,v)j  w# (U 2>  <  w,  v2  <  -  ^  /w*(u  )eXp(^  C)f(C,U2)dC>  . 

Now  for  0  <  U1  <  U2,  **(1),)  <  w*(U2)  and  f(C,U)  *  p(C)  -  p(w_)  + 

U2(?-w_)  so  f(Z,U2)  >  f(C,U^).  Therefore  if  (w,v)  €  8y  we  have 

°2 

w#(U.j)  <  w*(U2)  <  w  and 


v 


2 


< 


< 


fC(u2)e^lf  0£lt-u2)* 

lC<U,)« P<f  <>««-»,>*  • 


Hence  (w,v)  €  8  . 

^  a 


□ 


Lemma  4.9.  Assume  there  exists  U  >  0  so  that  there  is  a  solution  of  (3.8) 
connecting  the  saddle  (w_,0)  with  the  saddle  (w*(U),0)  in  the  v  >  0  half 
plane,  w  (U)  <  w+(U).  Then  there  exists  e  >  0  so  that  if  -e  <  u-U  <  0, 
Yy(w_,0)  has  (w+(U),0)  for  its  w-limit  set. 

Proof .  For  the  orbit  Yy(w*(U),0)  exiting  from  the  saddle  (w*(U),0)  in 
the  v  >  0  half  plane  set  y(w,U)  =  w* (5).  Then  y  satisfies 

A  y  ^  +  Uy  +  Dy2  +  p(w)  -  p(w_)  =  0  ,  (4.7) 

y(w*(U),U)  =*  0  .  (4.8) 
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From  (4.7)  we  have 


y2(w,u)  -  -  7  U  lV  _  exp<~  (C-w)  )y(C,U)« 
w*  (U) 

-  |  /W  exp(^  (C-w))f(C,U)<JC  . 

A  w*(U) 

Thus  for  w* (U)  <  w  <  w+(U)  we  have 

n2(w,o)  -  y2(w,u)  -  n2(w*(u),u)  =» 

-  I  /W  exp(~  (C-w))(un(C,U)  -  Uy(C,U))d£ 

A  w*(U) 

-  |  /W  exp(~  l?-w))(f(5,U)  -  f(C,U))dC 
A  w*(U) 

where  n(w,U)  satisfies  (4.3),  (4.4).  Thus  we  find 
h2(w,U)  -  y2(w,U)  -  n2(w*(U),U)  * 

-  |  (U-U)  exp(~  (C,w))n(C,U)d C 

A  w*(U>  A 

-  |  U  /W  exp(~  (C-w) ) (n(C,U)  -  y(C,U))dC 

w*(U) 

-  |  (U2-02)  /W  _  exp(^  (C-w)) (C-w  )dC  . 

w*(U) 

Now  if  n(z,U)  -  y(z,U)  <  0  for  some  z,  w*(U)  <  z  <  w+(U),  we  have 

Y*(w_,0)  entering  fl_  c  fl  and  we  can  use  Lemma  4.6  to  obtain  the  result. 
°  0  U  _  _ 

On  the  other  hand  if  n(C,U)  -  y(C,U)  >  0  on  [w*(U),  w+(U)]  then 

n2(w,U)  -  y2(w,U)  -  n2( w(U),U)  < 

-  j  (U-U)  /W  _  exp(~  (C-w)  )n(C,U)dC  +  const.  |U2-U2|  , 

w*(U) 

where  const,  depends  only  on  w*(U)  and  w+(U).  From  Lemma  4.4 
n(C,U-e)  >  n (C ,U)  >  0  on  (w_,w(U)).  Since  we  assume  n(C,U)  >  y(C,U)  on 
(w*(U),  w+(U)]  and  we  know  y(C,U)  >  0  on  (w*(U),  w+(u)]  we  must  have 
w  (U)  <  w(U).  Thus  (w  ,w(U))  must  contain  (w*(U),  w  (U)]  and  therefore 


n2(w,u)  -  y2(w,u)  -  n2(w(u),u)  < 


-  ^  (U-U)  Iw  _  exp(^  (C-w))n(C,U-e)dC 

A  w*(u)  A 

+  const. |U2-U2|  <  const. |U-U|  , 


(4.9) 


on  [w*(U),  w+(U) ] •  So  if  we  can  show  n2(w*(U),U)  <  const. |U-U|  we  will 

have  0  <  h2(w,U)  -  y2(w,U)  4  const. |U-U|  on  [w*(U),  w+(U)J.  Hence  for 

-e  <  u-U  <  0  and  e  >  o  sufficiently  small,  Y*(w  ,0)  must  enter  Q_  c  fl 

u  -  u 

and  Lemma  4.6  will  imply  the  result.  As  advertised  above  to  conclude  the 

2  — 

proof  we  need  to  estimate  n  (w*(U),U).  From  (4.5)  we  have 
n2(w*(u),u)  -n2(w*(u),u)  -  n2(w*(u),u)  - 


-  |  U  J**(U)exp(jp  (C,w))n(i;,0)dC  -  |  J"*(0)exp(®(C,w))f(C,U)dC 
+  7  u  J"*(U)exp{^  (C-w))n(C,U)dC  +  |  /w* (U ) eXp(^.  (c-w)  )f  (C,u)dC 


A  '  w_ 
2  rw* (U) 


-  If 

A  1  W 


2D 


exp(f=1  (C-w))(Uh(C,U)  -  On(C,U))dC 
A 


-  |  (u2-u2)  J"*(U)ex p(jp  <C-w))(C-wJdC 
”  A  ^w#(”)eXp(?  (C~w) )  [  (U-U)i)(C»U)  + 
U(n(C,U)  -  n(C,U))]dC 


a> 


-  f  (u2-u2)  /”*(u)exp(^  (C-w)Mf-wJdC  . 


2D 


A  w_  A 

Since  n(C,U-e)  >  n(C,U)  >  n(?,U)  >  0  on  (w  ,w*(U))  we  see 


n  2(w*(U),U)  4  -  |  (U-U)  /"MU)ex p(|^  (C-w))n(C,U-e)dC 


2  .  2  ~ 2  fW*(U)  2D  ,  , 

-  —  (U  -U  )  J  exp(—  (C“w) ) (C-w_)dC 


2  —  — 

and  hence  h  (w*(U),U)  4  const. |U-U|.  The  proof  is  now  complete. 


□ 
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n 


Lemma  4. 10.  Assume  there  exists  U1  >  0  so  that  the  orbit  Yn  (w  ,0)  of 
(3.8)  has  (w+(U.j),0)  as  its  w-limit  set.  Then  there  exists  e  >  0  so 
that  if  -e  <  U  -O.  <  0  then  Y~  (w  ,0)  has  (w+(U2)*0)  as  its  w-limit 
set. 

Proof.  Since  fly  <=  and  Y*  is  continuous  with  respect  to  U,  we  know 

Y*  (w  ,0)  will  enter  ft  for  U2  near  Uj.  How  use  Lemma  4.6. 

U2  U2 

□ 

Lemma  4.11.  Assume  there  exists  U  >  0  so  that  there  is  a  solution  of  (3.8) 
connecting  the  saddle  (w_,0)  with  the  saddle  (w*(U) ,0)  in  the  v  >  0  half 
plane,  w*(U)  <  w  (0).  Then  if  A  >  7  we  have  for  all  U,  0  <  U  <  O, 

T  O 

Yy(w_,0)  has  (w+(U),0)  as  its  w-limit  set. 

Proof.  From  Lemmas  4.9  and  4.10  we  know  the  conclusion  of  Lemma  4.11  will  be 

true  unless  there  is  a  decreasing  sequence  {u  },  U  >  0  and  U*  >  0  so  that 

n  n 

O  \U*  as  n  ♦  •  , 

n 

(w  (u  ),0)  is  in  the  u-limit  set  of  Y*  (w  ,0)  , 

+  n  U  — 

n 

(w+(U*),0)  is  not  in  the  w-limit  set  of  Y^(w_,0)  . 

From  Lemma  4.2  we  know  n(w,U*)  is  bounded  on  (w_,w(U*))  by  some  constant 

K(U*).  By  Lemma  4.4  we  know  n(w,U  )  is  bounded  by  K(U*)  on  (w  ,  w(U  )), 

n  —  n 

U*  <  U  <  U,.  Since  w(U  )  *  w,(U  )  we  have  n(w,U  )  bounded  by  K(U*)  on 
n  l  n  +  n  n 

(w_,  w+(°n^*  Prom  (4.5)  we  have  on  (w_,  W+(U„)) 

o  <  n2(w,u*)  -  n2(w,u  )  - 

n 

7  /w  exp (7^  (C-w)H(U  -U*)(n(C,U*)  -  n(C,u  )) 

A  w_  An  n 

-  (u  n(c,u*)  -  u*n(C,u  ))]dc 
n  n 

-  7  (U*2-U2)  /w  exp(^  (C-w))(C-w  )dC  . 

A  n  w  A  * 

By  Lemma  4.4  U  n(?,U*)  -  U*h(C,U  )  >  0  on  (w  ,  w,  (u  ))  and  since  we  have 
n  n  -  f  n 
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the  bound  K(U*)  on  n(C,U#)  and  n(?,U  )  (both  of  which  are  positive  on 

n 

(w_/  w+(Un))  1  we  see  0  <n2(w»u*)  ”  n2(w,Un)  <  const.  |Un-U*|  on 

(w  ,  w  (U  )).  Now  let  w  ♦  w  (U  ).  Since  n2(w  (U  ),U  )  »  0  we  have 
-  +  n  +  n  *  n  n 

0  <  n2(w  (U  ),U*)  <  const.  |U  -U*|  and  so  n2(w  (U  ),U*)  ♦  0  as  U  ^  U*. 

+  a  n  +11  n 

Also  since  w  (U  )  ♦  w  (U*)  as  U  'si  U*  it  follows  that  (w,(U*),0)  is  in 
+  n  +  n  T 

the  w-limit  set  of  Yy*(w_,0).  This  is  the  desired  contradiction  and  the 
lemma  is  proven. 

□ 


Lemma  4. 1 2.  Assume  A  >  0  and  let  w_  be  such  that  either 
(I)  J^*(0,exp(-~£  <  0,  Y  <  w_  <  a  , 

or 

(I* )  b  <  w_  <  Y  , 

and 

(XI)  assume  there  exists  UM  >  0  so  that 

w*  (U  ) 

/„  M  exp(~  C)f(C,UM)dC  "  0  • 

Then  there  exists  a  U,  0  <  U  <  u  so  that  (3.8)  possesses  a  solution  in 

n 

the  v  >  0  half  plane  connecting  the  saddles  (w_,0)  and  (w*(U),0). 
Proof.  From  (4.5)  we  know 

n 2(w,u)  -  -  ?/w  exp(~  (c-w))n(5/u)dc 
A  W  A 

-  7  exp(-f  (?-w))f(C,U)d? 

AW  A 

where  n  satisfies  (4.3)/  (4.4). 

For  U  *  0, 

n2(w,0)  =  -  |  f”  exp(^  (C-w))f(C,0)dC  . 
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If  (I')  holds  f(C,0)  *  p(C)  -  p(w_)  <  0  and  n2(w,0)  >  0.  If  (I)  holds 

2 

(H1)-(H6)  imply  again  n  (w,0)  >  0.  So  the  orbit  exiting  (w_,0)  for  (3.8) 
in  the  v  >  0  half  plane  never  crosses  the  w  axis.  On  the  other  hand 

2U  w*  (U  )  __ 

n2(w*(UM),UM)  =  -  -j1  fw  M  exp(—  (C-w)  )n(C,UM)dC 
w*  (U  ) 

-  7  /  M  exp(-~  <C-w))f(C,Uj« 

AW  A  H 

»M  W*  ( U  ) 

=  -  -~/  exp(^  (C-w))n(C,UM)dC 

AW  A  n 

which  implies  n(C,U  )  =  0  for  some  z,  w  <  z  <  w*(U  ).  By  continuity  with 
respect  to  U  we  see  that  for  some  U,  0  <  U  <  UM,  n (w*(U) ,U)  *  0.  The  proof 
is  complete. 

□ 

Proof  of  Theorem  4.1.  Part  (i)  follows  from  Lemma  4.12.  Part  (ii)  is  a 
consequence  of  Lemmas  4.10/  4.11  and  4.12. 

□ 

We  now  give  a  mechanistic  interpretation  of  Theorem  4.1.  Let  us  assume 
A  >  and  w_  is  such  that  (I)  or  (I*)  and  (II)  hold.  If  we  look  on  the 
Huginot  curve  (Figure  9)  we  see  for  w  (U  )  >  w  (U),  0  <  U  <  U,  (4.2)  with 

T  &  T  A 

U  *  U.  gives  the  solution  of  the  Riemann  problem  (1.6)  with  Cauchy  data 

Cl 

(w  =  w  w  -  W+(U  ) 

X  <  0,  X  >  0  .  (4.10) 

u  *  u,  +  U  (w, (U  )-w  )  u  *  u. 

+  a  +  a  -  + 

This  solution  is  admissible  according  to  the  viscosity  -  capillarity 
criterion  (A,D),  A  >  —.  A  representation  of  the  solution  X  -  t  plane  is 
given  in  Figure  10. 
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(U  )  -  w  ) 

cl 


w  =  w  (U  ) 

t  d 


Figure  10:  Solution  of  (1.6),  (4.10) 

As  we  compress  the  fluid  so  that  w*(U)  <  w+ ( )  <  w+(U)  for  some  Ub» 
Ojj  >  0  >  0,  the  adsdssibillty  of  the  above  solution  is  no  longer  guaranteed 
by  Theorem  4.1.  For  example  we  see  the  Riemann  problem  for  Cauchy  data 


X  <  0, 


u  -  i  +  U,  (w(U.)  -  w*(U) ) 
+  D  +  D 

+  U(w*(U)  -  w  ) 


w  ■  »+(Ub) 


u  *  u. 


X  >  0 


(4.11) 


a  solution  (see  Figure  11) 


u  -  u+  +  Ub(w+(Ufa)  -  w*(U) )  +  U(w*(U)  -  w_)  X  <  Ut  , 


w  -  w*(U) 


u  »  u  +  U(w  (U.  )  -  w*  ( U ) ) 
T  T  D 


Ut  <  X  <  Ut  , 


(4.12) 


*  -  VV' 


u  «  u. 


Ut  <  X  , 


where 


J 

*  •  _ 

I 

I 


u. 


■p(w+(Ub)  )+p(w_) 
w+(Ufa)-w_ 


I 


U 


u 


J -p(w*(U))+p(w+(Ub)) 
w* (U)-w(U) 

+  D 

^/-p(w*(U))+p(wJ 


w*(U)-w 


t 


From  Theorem  4.1  and  Theorem  3.3  (i)  we  see  the  solution  is  admissible 


according  to  the  viscosity-capillarity  condition  (A,D) ,  A  >  —  in  that  both 

O 

X  *  Ut  and  X  =  Ut  are  admissible  singular  surfaces.  Such  a  solution 
indicates  that  as  the  fluid  is  compressed  for  w  (U  )  >  w  (U)  a  single 

’  3L  • 

propagating  phase  boundary,  separating  liquid  and  vapor  phases,  and  traveling 
with  speed  Ua  >  0  is  possible.  However  upon  further  compression 
w+(U^)  <  w+(U),  this  interface  bifurcates  into  a  classical  shock  in  the  vapor 
phase  moving  out  from  the  propagating  interface.  The  shock  travels  with  speed 

<v  —  — 

U,  the  liquid-vapor  interface  now  travels  with  fixed  speed  U.  Thus  U 
represents  the  maximum  speed  a  liquid-vapor  interface  may  travel  in 


On  the  other  hand  if  we  reverse  the  process  and  start  with  initial  data 


(4.11)  and  raise  w^U^)  through  expansion  the  shock  wave  and  interfacial 
boundary  of  (4.12)  will  coalesce  into  an  interfaciil  boundary  traveling  with 
speed  Ufl. 

Example  1 2.  Let  w+  be  given  in  (0,“)  and  let  w_(U)  be  a  solution  of  the 
Rankine-Huginot  relation  with  U  <  0,  namely 


U 


^ -p(w+)+p(w_(U) ) 


w+-w_(U) 


In  this  case  an  argument  similar  to  that  given  for  Lemma  4.12  yields  the 
following  result. 

Theorem  4.13.  Assume  A  >  0  and  let  w+  be  such  that  either 
w 

(I)  /w+(0)exp(j—  C)f(C,0)dC  >0  0  <  w+  <  6  , 

or 

(I*)  6  <  w+  , 

and 

(IX)  assume  there  exists  UM  <  0  so  that 

w 

C(u  )exp(5Tc)f(?'0)dc  “  0  * 

M 

Then  there  is  a  unique  U,  U  <  U  <  0,  so  that  (3.8)  possesses  a  solution  in 

N 

the  v  <  0  half-plane  connecting  (w_(U),0)  to  (w+,0). 

Pictorially  this  result  is  represented  in  Figures  13  and  14.  In  this 
case  the  singular  surface  X  *  Ut  separating  liquid  and  vapor  phases  is  a 
rarefaction  shock  which  is  admissible  according  to  the  viscosity-capillarity 
condition  (A,D).  This  example  stands  in  contrast  to  Example  9  where  the 
interfacial  rarefaction  shock  was  not  (modulo  exceptional  circustances ) 


admissible  according  to  the  viscosity-capillarity  condition  (A,D) 


Figure  14:  Rarefaction  shock  waves 


<U) 
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If  the  fluid  is  compressed  or  expanded  at  changing  to  value  of  w_ 

the  pure  rarefaction  will  disappear.  Expansion  will  set  new  Cauchy  data 
if  w  (U)  <  w  <  a, 

—  A 

w»w  X  <  0  ,  w  *  w^  X>0 

a  + 

_  (4.13) 

u  *  u  u  ■  u  +  U(w  (U)-w  )  +  U(w  -w  (U)) 

•  *■  ™  fit  T  * 


then  (1.5)  possesses  a  solution 


On  the  other  hand  if  we  compress  the  fluid  at  + 


we  can  envisage 


Cauchy  data 


w  «  w.  X  <  0,  w  =  wJ  X>0 
D  + 

u  *  u_  u  *  u+ 


(4.15 


where  b  >  w^  <  w_(U).  In  this  case  we  again  can  connect  the  vapor  state 
w  *  w+,  u  ■  u+  to  a  liquid  phase  w  *  w_(U),  u  =  u+  -  U(w+-w_) (U) ) .  Then 
this  liquid  phase  can  be  connected  by  a  classical  rarefaction  wave  in  the 
liquid  phase  to  w  =  wb,  u  **  u_  for  w^,  u_  sufficiently  near 
w_(U),  u+  -  U(w+-w^(U))  (see  Lax  [16]).  Pictorically  the  result  is  shown  in 
Fiqure  16. 


Figure  16 

Example  13.  In  this  example  we  see  that  "metastable"  states  are  dynamically 

admissible  as  perturbations  of  f  an  equilibrium  configuration.  Assume  u. 

u+  =*  0,  U  *  0,  and  w_  and  w+ 

w 

/w+exp(~  i  )  (p(C  )  -  p(w_))d?  »  0  . 


Hence  according  to  Example  10  this  singular  surface  is  admissible  according  to 


I 

1 

j 

I 


the  viscosity-capillarity  condition  (A,D) .  But  from  Theorem  3.3  (i)  the 
singular  surfaces 

T  :  X  -  U  t  separating  w  *  w  ,  u  *  0  and 

w  “  w  ,  u  -  -U  (w, -w  )  with  B  <  w  <  w,  and 
a  a  +  a  a  + 

/-p(w  )+p(w  ) 

u  -  +  /  — - - f 


r  ^  t  X  *  U^t  separating  w  *  w_ ,  u  -  0  and 

w  “  w,  ,  u  •  -U.  (w  -v.  )  with  w  <  w.  <  a  and 
D  D-D  “  D 


^/-p(wfa)+p(w_) 

V*- 


are  admissible  according  to  the  viscosity-capillarity  criterion  (A,D),  A  >  0. 
Thus  the  Cauchy  problem  (1.6)  with  initial  data 

w  *  w  w  ■  w 

"  A 

x  <  0  ,  X  >  0  (4.16) 

u  *  0  u  ■  -0  (w  -w  ) 

A  »  A 

and 


w  =•  w, 

D 

u  *  -U.(w_-w  ) 


X  <  0 


u  ■  0 


X  >  0 


possesses  solutions 


(4.17) 


w  *  w_  w  *  w+ 

X  <  0, 


u  *  0 


u  =*  0 


w 

0  <  X  <  U  t, 

a 

u 


X  >  U  t 
a 


-U  (w  -w  ) 

A  *  A 


(4.18) 
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I 


,1 


1 

» 

i 

l 


b  b  -  '  . 

u  -  -0.  (w  -v.  )  u  »  0  u  *  0 

JO  —  D 

By  Theorem  3.3  (i)  and  Example  10  these  solutions  are  admissible  according 
the  viscosity-capillarity  condition  (A,D) ,  A  >  0.  (See  Figure  17.) 


In  both  examples  we  see  solutions  exist  in  the  "metastable"  regions 


(w_,a)  and  (0,w+)  though  only  for  a  period  of  finite  duration  determined 
by  U.  and  Uh  respectively. 


Remarks .  The  results  presented  In  this  section  are  consistent  with  those 
given  by  Zel'dovich  and  Raizer  [17,  p.  750-762].  However  unlike  the  purely 
physical  arguments  given  there  our  results  are  based  on  a  rigorous  analysis  of 
the  "shock  structure  problem".  Moreover  we  note  that  the  result  of  Theorem 
4*1  (i)  showing  the  theoretical  existence  of  a  shock  wave  (admissible  in  the 
sense  of  the  capillarity-viscosity  criterion  (A,D),  A  >  0)  for  the  complete 
transition  from  superheated  vapor  to  liquid  is  consistent  with  the 
experimental  results  of  Dettleff,  Thompson,  Meier,  and  Speckman  [18];  see  [5] 
for  a  complete  non- isothermal  discussion  of  this  result. 
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